Survival probability of a diffusing particle in the presence of Poisson-distributed 

mobile traps 
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The problem of a diffusing particle moving among diffusing traps is analyzed in general space 
dimension d. We consider the case where the traps are initially randomly distributed in space, 
with uniform density p, and derive upper and lower bounds for the probability Q{t) (averaged 
over all particle and trap trajectories) that the particle survives up to time t. We show that, for 
1 < d < 2, the bounds converge asymptotically to give Q{t) ~ exp( — A^f'*''^) for 1 < d < 2, 
where Xd = (2/7rci) sin(7rci/2)(47r_D)'*''^p and D is the diffusion constant of the traps, and that 
Q{t) ~ exp{—4npDt/ \nt) for d — 2. For d > 2 bounds can still be derived, but they no longer 
converge for large t. For 1 < d < 2, these asymptotic form are independent of the diffusion constant 
of the particle. The results are compared with simulation results obtained using a new algorithm 
[V. Mehra and P. Grassberger, Phys. Rev. E 65, 050101 (2002)] which is described in detail. 
Deviations from the predicted asymptotic forms are found to be large even for very small values 
of Q{t), indicating slowly decaying corrections whose form is consistent with the bounds. We also 
present results in d = 1 for the case where the trap densities on either side of the particle are 
different. For this case we can still obtain exact bounds but they no longer converge. 

PACS numbers: 05.40.-a, 02.50.Ey, 82.20.-w 



I. INTRODUCTION 

Reaction-diffusion processes represent a large and im- 
portant class of systems with nonequilibrium dynamics. 
From a fundamental physical viewpoint, the interest in 
these systems lies in the fact that the concentration of 
reactants is governed, in general, by irreversible reaction 
events that depend on the spatial distribution of par- 
ticles rather than through equilibrium fluctuations con- 
trolled by a chemical potential. Such model systems have 
a range of applications, most notably to chemical kinetics 
[ 1, 01 but also to interfacial growth domain coarsening 
[ 4 , g and aggregation . 

The most intensively studied reactions are single- 
species annihilation {A + A ^ ^) and coalescence {A + 
A—^ A) as well as two-species annihilation {A+B 0) — 
see, e.g., @, |[ |§1 for reviews. In this paper we focus on 
the two-species problem. It is known to exhibit two differ- 
ent classes of long-time behavior depending on whether 
the initial concentrations of A and B particles are equal 
or not. (As an aside, we note that a similar dependence 
on the initial condition also holds for the A -I- A ^ re- 
action when the reactant motion is deterministic rather 
than diffusive The reason for this is that when 

the initial densities of A and B particles are the same, 
they remain so for all times, whereas if, say, the initial 
density of A particles pa(0) is less than that of the B par- 
ticles pb{0), the ratio pa{0)/pb{0) ^ as i ^ oo and at 
late times one has a few, isolated A particles diffusing in 
a background of B particles. 

The case of equal initial densities is well understood, 
and results similar to those for the A + A ^ with diffu- 
sive particle motion have been obtained [|l2|, In low 



dimensions, here d < 4, fluctuation effects are important 
and one finds a density decay pAit) = Psit) ~ t"**/^ 
in this diffusion-limited regime. Above the critical di- 
mension d > = 4 one finds that the mean-field result 
PA,B{t) ^ 1/i applies. This result also holds for the 
A + A ^ % process above its critical dimension dc — l. 

By contrast, the density decay forms for the A+B — > 
process when the initial densities pA{t) and psit) are not 
equal are less well understood. In fact, since the exposi- 
tion of the process as a model of monopole-antimonopole 
annihilation in the early universe nearly twenty years ago 
|l2| , only a few results are known exactly. Most notably, 
Bramson and Lebowitz proved rigorously that, at 
large times, the density of the minority species (which 
we will take to be the A particles) behaves as 



PA{t) 
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d^2 
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(1) 



revealing d = 2 to be critical in this case. To the best of 
our knowledge, no predictions for the constants A^ were 
given until recently [^5j. Furthermore, there has been 
no convincing numerical verification of the predicted de- 
cay even in one-dimension, despite the development of 
sophisticated simulation techniques jl^ that allow the 
probing of extremely small densities that emerge at large 
times. In this paper, we expand on the bounding argu- 
ments reported in that give rise to precise values of 
Xd for d < 2. We also present a detailed description of 
the simulation algorithm introduced in and extend 
it to test our bounding arguments and understand the 
approach to the asymptopia described by Eq. (^. 

As noted above, the late-time regime is characterized 
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by a few isolated A particles diffusing in a sea of B parti- 
cles. Thus it is appropriate to consider the extreme case 
of a single A particle in a sea of B particles that has a uni- 
form (Poisson) density. In this case, the quantity pA{t) 
is just the survival probability of the A particle. Further- 
more, if the diffusion constants of the A and B particles 
are the same, one can also view pA{t) as the fraction of 
particles that have not met any other particles. Thus the 
reaction A + B in the limit of a low density of A par- 
ticles has been discussed under the guises of uninfected 
walkers in which random walkers infect each other 
on contact, diffusion in the presence of traps ]l6| , [T8| in 
which the B particles are considered as traps for the A 
particles, and predator-prey models ]l9| in which one asks 
for the survival of a prey (the A particle) being 'chased' 
by diffusing predators (the B particles). To avoid con- 
fusion, we shall adopt only the trapping terminology in 
our discussion. 



In this work, we show how the survival probability of a 
diffusing particle in the presence of mobile traps can be 
understood in terms of the target annihilation problem 
pO| , [2l| , (or first passage problem [^) where one asks 
for the probability that none of the traps has entered a 
particular region (target) in the d-dimensional space. In 
turn, the asymptotics of the target annihilation problem 
are intimately related to the recurrence or transience of 
diffusion in various dimensions. A process is said to be 
recurrent if the probability of returning to the initial con- 
figuration is unity: in the context of diffusion, this implies 
that with probability one a walker will visit a particular 
point in space infinitely often. It is well known (see, e.g, 
p3| , |2^ ) that diffusion is recurrent in dimensions d < 2, 
whereas in more than two dimensions it is transient (i.e. 
the return probability is less than one). It is precisely 
this property of diffusion that gives rise to the critical 
dimension of two for the trapping reaction and hence the 
asymptotic results (^ for the A + i? ^ process. 

The principal result of the paper is the determination 
of the constants in Eq. (|l]) for d < 2, and the deriva- 
tion of upper and lower bounds for d — 3. A striking 
feature of the results is that, for d < 2, the value of is 
independent of the diffusion constant of the A particle. 

We begin in the next section of this paper by defin- 
ing the trapping reaction model. Then, in section [II 
we present in detail our analysis of the one-dimensional 
case, testing our predictions in section ^ where we dis- 
cuss how the model may be simulated efficiently. In sec- 
tion we show how the method used to treat the one- 
dimensional case can be extended to general dimensions 
d > 1. Only when the underlying diffusion process is re- 
current (i.e. for d < 2) do our upper and lower bounds 
converge asymptotically to give exact predictions for A^. 
Finally, in section we present a discussion and sum- 
mary of the results. 



II. DEFINITION OF THE MODEL 

The trapping reaction model we consider is defined 
as follows. At time t = a particle is placed at the 
origin of a d-dimensional coordinate system. Surround- 
ing this particle is a uniform sea of traps whose initial 
positions Xi are chosen independently. This initial con- 
dition ensures that the distribution of traps is Poisson, 
i.e. the probability that a volume V contains N traps is 
[{pV)^ /Nl] exp{—pV) in which p is the mean number of 
traps per unit volume. 

The dynamics of the particle and traps can be ex- 
pressed using the Langevin equation 



(2) 



in which the subscript i = denotes the particle, i > 1 
one of the traps and the superscript a indicates a com- 
ponent of the position vector Xi. The noise rif{t) is a 
Gaussian white noise with zero mean and correlator 



ivTiMit')) ^2DAjS^pS{t-t') 



(3) 



We take all the traps to have a diffusion constant D and 
the particle to have a diffusion constant D' . Hence Dq — 
D' and Di = D ioi i > 0. The quantity of interest in 
this model is the probability Q{t), averaged over all initial 
conditions and realizations of the random walks, that the 
particle has not yet met any of the diffusing traps. 



III. ANALYTICAL RESULTS IN ONE 
DIMENSION 

For clarity, we restrict ourselves initially to the case 
d = 1. Later, in section ^ we will explain how the ar- 
guments presented in detail here can be generalized to 
higher dimensions. We begin with a description of the 
target annihilation problem before moving on to discuss 
how it applies to the more general problem of a particle's 
survival in a sea of diffusing traps. The target annihila- 
tion problem can be solved exactly for any d ^ . 
The asymptotic form of the solution, and the leading 
corrections to it (for d > 1), play a central role in our 
bounding arguments. To establish the notation and to 
make our presentation self-contained, we present in this 
paper a brief derivation of the main results as a prelude 
to deriving the bounds. 



A. The target annihilation problem 

Consider a one-dimensional line containing a target of 
length 21 centered on the origin (i.e. lying between x — ^l 
and X = I). We wish to calculate the probability Qrit) 
that none of the diffusing traps initially placed outside 
this region has hit the target by a time t. This quantity 
can be calculated if one knows the probability Qi{t\y) 
that a trap initially at position y has not yet entered the 
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target region. Since the target is static and each trap 
executes independent diffusion, we can simply multiply 
the probabilities for each individual trap together and 
average over all possible initial positions to find Qxit). 

Let us consider then a trap that has its initial position 
to the right of the target, i.e. y > I. The probability 
Qi(t\y) that the trap has not reached the target satisfies 
the backward Fokker-Planck equation 



dQi{t\y) Qi(%) 



dt 



9y2 



(4) 



with the boundary conditions Qi{t\l) = 0, Qi(0|y) = 1 if 
y > I and Qi(t\oo) = 1. These express the facts that the 
probability that the target has been reached if the trap 
started at y = Z is one, that it is reached in zero time 
from y > I is zero and that it is reached from infinity in 
a finite time is zero respectively. The solution to (Q) that 
satisfies these boundary conditions is 



Qi(%) = crf 



(5) 



in which erf (x) is the error function. 

Instead of a single trap to the right of the target, 
consider N independently diffusing traps, each initially 
placed at random in the interval yi G [1,1 + L]. Then, the 
probability that none of the traps has reached the target 
by time t is 



IN 



N -I A+L 



dyi erf 



(6) 



It is convenient now to rewrite the error function in terms 
of the complementary error function, erf (a;) = 1 — erfc(a;). 
Then one has 



QNit) 



l+L 



dy erfc 



y 



I 



N 



(7) 



Since we wish to consider an infinite sea of traps, we take 
N — pL and then the limit L — > oo holding p, the density 
of traps, fixed. This yields 



Qooit) = lim 

L — >oo 



exp 



1 

1 

L 



dy crfc ^ 



pL 



2pVDi 



(8) 



This gives the probability that no traps initially posi- 
tioned on one side of the target have reached the target 
by time t. Since we have in mind a target surrounded 
on both sides by traps, and that the motion on each side 
is independent, we obtain the probability that the target 
has not been annihilated by a trap by squaring (0) . That 



IS, 



Qrit) = exp 



ApVDt 



(9) 



Note that the size of the one-dimensional target I does 
not appear in this exact expression for its survival prob- 
ability. Later, in section ^ we will find that at suitably 
large times, the size of the target is unimportant for all 
d < 2 (where diffusion is recurrent). 



B. Bounding argument for a diffusing particle in 
the presence of mobile traps 

We now discuss how to construct upper and lower 
bounds on the particle's survival probability Q{t) us- 
ing the result for the target annihilation problem (S) in 
one dimension. We claim that, on average, a particle 
surrounded by a uniform, isotropic distribution of traps 
survives longer if it is stationary than if it is allowed to 
diffuse. We are currently unable to prove this statement 
rigorously, although it is supported by intuition and nu- 
merical data (see section ^ below) . We also note that 
when we say "on average" we mean "after averaging over 
all possible initial trap positions and trajectories of both 
particle and traps" . 

If this claim is accepted, we obtain an upper bound 
Qu{t) on the particle's survival probability from (^ by 
noting that requiring the particle to remain stationary is 
equivalent to having a target region of size 1 = 0. Then 
we immediately have that 



Q{t) < Quit) = exp 



ApVDt 



(10) 



To derive a lower bound on the survival probability 
Q{t) we introduce a notional box of size / centered on 
the origin. If we ask for the particle to remain inside this 
box until a time t, and for all the traps to remain out- 
side it, the traps and particle may never meet and hence 
the particle survives until time t. There are, of course, 
other trajectories for which the particle survives, and so 
those just described form a subset of all possible surviv- 
ing trajectories — see Fig. ^ Hence the probability that 
the particle remains within the box and traps outside is 
a lower bound QL{t) on Q{t). 

There are three independent contributions to this 
bound: (i) the probability that there are initially no traps 
in the box of size /; (ii) the probability that no traps enter 
the box up to time t; and (iii) the probability that the 
particle has not left the box up to time t. The first two 
contributions are easily obtained. From the definition of 
the Poisson distribution, we have that the probability the 
box initially contains no traps is exp(— pZ). Secondly, the 
probability that no traps enter the box is independent of 
the box size and is given by (^. The third contribution, 
the probability that the particle remains inside the box, 
is obtained as follows. 

Since the system is translationally invariant, we can 
just as easily consider a particle initially sandwiched be- 
tween absorbing boundaries at a; = and x = I. The 
probability Qp{t\y) that the particle starting at y = 1/2 
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(i) 




FIG. 1: Two walker trajectories (space-time plots, with t — 
at the bottom) for which the particle (unfilled) survives 
contact with a trap. Trajectories of type (i) have the property 
that the particle remains inside a notional box, and the traps 
outside. This forms a subset of the entire class of surviving 
trajectories, which includes paths of type (ii) in which the 
particle leaves the box and the traps enter but nevertheless 
no particle-trap contact occurs. 



once the time t is sufficiently large. Note that this pro- 
vides a bound for a particular box size I. Since the box 
is an artificial construct, we can choose its size so that 
the lower bound is maximized at a particular (predeter- 
mined) time t* . One finds that the corresponding box 
size is I* = {2T:'^D't* / pf/^. Using this box size in (|l|) 
we find the largest lower bound is given by 



(16) 



QhKt) = - cxp -= 3 

TT \ Jit 





l/3\ 
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Combining this lower bound with the upper bound 
Quit) of Eq. we find 



4 ^ Ing(t) 

- (p^Dty/^ 



< 



-^[2) {p^Dty/^ • ^^^^ 



This implies that the constant Ai in the expressions of 
Bramson and Lebowitz (0) is precisely determined as 



has not crossed the absorbing boundaries satisfies a back- 
ward Fokker-Planck equation 



dQp{t\y) _ ^,d^Qp{t\y) 



dt 



= D 



(11) 



subject to the absorbing boundary conditions Qp{t\Q) = 
Qp{t\l) — and the initial condition Qp{t\y) — 1 for 
Q < y < I. The general solution to this equation that 
satisfies the absorbing boundary conditions is obtained 
by separating the time and space variables in the usual 
way. One obtains the Fourier sine series 



Qp{i\y)^Y.^^'''^^\ 12 j'"\~rj ■ ^^^^ 



The expansion coefficients are fixed through the initial 
condition. Using the orthogonality of the sine functions 
one finds 



— k odd 
k even 



(13) 



For the purposes of the present calculation, we need 
consider only the long-time form oiQp(t\y) for a particle 
that starts at y = //2. Thus we keep only the longest- 
lived (fc = 1) mode in the expansion (ll^) to find 



Qp(t\ll2) - -exp 

TT 



(14) 



Ai 



lim 

t — ^00 



InQ(t) 



(18) 



Note that this constant depends only on the density and 
diffusion constant of the traps, and is independent of the 
diffusion constant of the particle. 



C. Extensions to the basic trapping reaction model 

It is straightforward to incorporate two generalizations 
of the one-dimensional trapping model defined in section 
^ into the bounding arguments discussed above. The 
first of these is to allow the traps to the left and right of 
the origin at time to have different densities. We denote 
the larger (respectively, smaller) of these densities as p+ 
{p^ ) and their average as p — ^ (p_|_ + p_ ) . Additionally 
we shall place n particles at the origin at time and 
study the probability that all survive until a time t. 

To obtain an upper bound on the survival probability, 
we note that the survival probability of the particles can 
only increase (or remain constant) as either or /3_ is 
decreased. Hence the survival probability for the case 
of unequal densities is bounded from above by the case 
where the density of traps is on both sides equal to p_. 
For the case of a single diffusing particle, we argued above 
that an upper bound on its survival probability is found 
by setting its diffusion constant D' to 0. Clearly, if Z?' = 
the number of particles at the origin is irrelevant, and so 
an upper bound on Q{t) is given by Eq. ( p^ with p — p-, 
i.e., 



Including this along with the contributions to the lower 
bound QL{i) on the diffusing particle's survival proba- 
bility Q(t) discussed above, we have 



Q{t) > Qdt) 



■ exp 



pi 



ipVDt 



(15) 



Q{t) < Q[/(t) =exp 



Ap-VDi 



(19) 



To obtain a lower bound on the particles' survival prob- 
ability we once again introduce a notional box, inside 
which all the particles must remain and no traps may 
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enter until time t. This time, however, we respect the 
asymmetry of the problem by allowing the box to extend 
a distance I- into the low-density region of traps and 
Z+ into the high-density region. We will again seek to 



J 



maximize the lower bound by varying Z_ and Z+. 

A lower bound Ql(^) is obtained using an argument 
analogous to that leading to Eq. Considering once 

again late times, we find 



Qhit) (X exp 



mr'^D't 



{p-i- +P+1+) 



(20) 



The number of particles n enters into this expression 
through the fact that the probability for all of the n 
particles to remain inside the box of size I = I- -\- 1+ is 
simply the v}^ power of the corresponding probability for 
a single particle. 

The maximal lower bound for a prescribed time t* is 
obtained from (20) by setting to zero (thus discounting 
particle trajectories that enter the high-density region) 
and putting /I = {2m:'^ D't* / p^f/^ . Then 



'mr'^p'iD't' 


l/3\ 
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Q{t) > Qdt) oc exp i-^P—Ei - 3 



Along with the upper bound ([l9| ) we find that 



(21) 



< -- 



InQ 



{plDtyi^ 



< 



pn p- 
4D 



1/3 



(22) 



Note that, except for the case where p — p^ (which im- 
plies p- ~ p+) these two bounds do not converge and 
so we cannot make a precise prediction for A when the 
trap densities are unequal. For the case p- = p+, how- 
ever, the bounds converge to 4/^/??, independent of the 
number of particles n. 



IV. SIMULATION ALGORITHM AND RESULTS 

A sophisticated algorithm for simulating the trapping 
reaction in discrete space and time and with a Poisson 
distribution of traps was recently introduced The 
beauty of the algorithm is that it admits (numerically) 
exact calculation of the survival probability for an arbi- 
trarily long, but fixed, trajectory of the particle. As will 
be discussed below, the algorithm takes into account all 
possible paths of the traps, as long as their initial dis- 
tribution is Poisson. In order to obtain an estimate of 
the particle survival probability, it is necessary to iterate 
the algorithm over a sequence of particle paths. We now 
discuss this algorithm in detail. 



A. An efficient simulation algoritlim 

In order to simulate the trapping reaction model in one 
dimension, we construct a discretized version in which 
each walker follows a path x{t) that has x{t + l) — x{t) = 
±1. Since all hops to the left or right occur in parallel, we 
must ensure that the initial coordinates of all the walkers 
are even integers so that no two walkers are able to hop 
over each other. 

As a starting point in understanding the simulation 
algorithm, consider a system comprising the particle, 
whose trajectory xo(i) is predetermined, and a single 
trap, whose trajectory xi{t) is stochastic given some ini- 
tial condition a;i(0) = t/i. The probability Pi(x,t) of 
finding the trap at site x after time t, given that it has 
not absorbed the particle, satisfies the equation 

Pi(x,t+1) = i[Pi(a:- l,t) + Pi(x + l,t)] (23) 

subject to the initial condition P{x,0) = Sx.yi and the 
moving absorbing boundary condition P{xo(t),t) — 0. 
Note that (|2^) is the discrete analogue of the diffusion 
(Fokkcr-Planck) equation 



dPiix,t) ^d^Pi{x,t) 



dt 



D- 



dx'^ 



(24) 



By Taylor expanding ( p3|) we find the diffusion constant 
of both particle and trap to be Z) = D' = i. 

The solution of the diffusion equation with an arbi- 
trary moving absorbing boundary at xo(i) is not known 
analytically. One can obtain it numerically, however, by 
iterating the following two steps over t' = 1, 2, . . . , t. 

1. Construct the probability distribution of the 
trap's position using the equation Pi{x,t') = 
i[Pi(x- l,t'- l) + Pi(x + l,t'- 1)]. 

2. Enforce the absorbing boundary condition by sub- 
sequently setting Pi{xo{t'),t') = 0. 

In the simulation, we wish to consider not just a single 
trap, but a Poisson distribution of traps. This can be 
achieved as follows. Let Pn{x,t) be the probability that 
there are n traps on lattice site x at time t. We shall 
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assume that this distribution is Poisson, i.e., 

Pnix, t) = M^lllT^ exp[-c(x, t)] (25) 
n\ 

in which c{x, t) is the mean number of traps at site x and 
time t. 

Now, if each trap can hop with equal probabihty to the 
left or right in one time step, we have 

n 

P„(x,t+1)= Y.W+(x-l,t)W-_„,(x + l,t) (26) 

m=0 

in which W^{x, t) is the probability that m particles hop 
from site x at time i to a; ± 1 at time < + 1. This quantity 
is given by 

^^(-^*)=f:^^exphc(.,t)]Q)l. (27) 

Insertion of this expression into (^6|) and a little algebra 
reveals that 

P„(x,t+l) = te^exp[-c(a;,i)] (28) 

in which 

-c[x,t) = '-^^^^2^±^^^±^. (29) 

That is, if the distribution of traps at time t is Poisson 
the distribution of traps at time t + 1 is also Poisson, 
with the mean occupation number at each site obeying 
the discrete diffusion equation 

c{x, t + 1) i [c{x - 1, i) + c(x + 1, t)] . (30) 

As with the case of the single trap described above, we 
wish to determine the probability distribution of traps 
given that the particle following the predetermined path 
X[){t) has not been absorbed until a time t. We must 
therefore have at each time step Pn{xo{t), t) — Sn.o which 
can be achieved by enforcing the boundary condition 
c{xo{t),t) = 0. Thus we can evolve the mean occupa- 
tion numbers for the Poisson distributed sea of traps in 
exactly the same way as for the single-trap distribution 
function described above (albeit with a different initial 
condition, to be described shortly). 

In the simulations, we wish to calculate the probability 
that the particle has survived until time t. To obtain an 
expression for this, consider a particular distribution of 
traps described by the function c(x,t). The probability 
that site a;o contains no traps is just exp(— c[a;o, t]) and so 
Q{t + 1) = Q{t) exp(— c[a;o, t]) where the value of c{xo, t) 
used is that obtained after the diffusion step, but before 
enforcing the boundary condition c(x{j,t) = 0. 

We now give a step-by-step explanation of the algo- 
rithm for calculating the particle survival probability for 



a predetermined particle path X{){t). One begins by set- 
ting up the trap concentration as follows: 

( 2pL X < xa{0) 
cix,to) = < a; = a;o(0) (31) 
[2pR x> xo{0) 

in which and pR are the equivalent continuum densi- 
ties t o the left and right of the particle, as used in section 
III C| . The factor of 2 emerges because that is the effec- 
tive lattice spacing in the discrete model. We also set 
(5(0) = 1 (i.e., we assume there are no traps at the origin 
to begin with). Then, for each time t' — 1, 2, ... ,t we 
perform the following steps: 

1. The trap concentration variables are evolved using 
c{x, t') = i [c{x - 1, t' - 1) -f c{x + l,t' - 1)]. 

2. The cumulative particle survival probability is cal- 
culated using Q{t') = Q{t' - l)exp[-c{xoit'),t')]. 

3. The boundary condition is enforced by setting 
cixoit'),t')=0. 

Note that this algorithm can be run for paths of arbitrary 
length and that, at a particular time the trap density 
at positions x < a;o(0) — t' and x > a;o(0) + 1' is uniform. 
Hence at each time step, one need deal only with t' + 1 
concentration variables to simulate the infinite system. 

Using the above algorithm, one obtains the survival 
probability for a particle following a particular path 
xoit). To reach an estimate of the particle survival prob- 
ability averaged over all paths, it is most efficient to per- 
form Monte Carlo sampling. That is, one generates a 
binomial random walk by choosing the particle displace- 
ment xo{t') — xo(t' — 1) = {—1, 1} with equal probability. 
Then, one estimates the mean particle survival probabil- 
ity as 

1 ^ 

QW«^E^?^'^W (32) 

fe=i 

in which Q'^'^^ is the value of the survival probability for 
the k^^ random walk. One can, of course, estimate other 
quantities, such as the mean and variance of the parti- 
cle's displacement. Also, if one is interested only in the 
short-time behavior, one can obtain the particle survival 
probability for each possible path. We should also note 
that the one-dimensional algorithm described here gen- 
eralizes straightforwardly to higher (integer) dimensions. 

B. Numerical results 

We first investigate the entire set of short particle paths 
in order to get a feel for those that give rise to the great- 
est probability of survival. For each time t < 28 we found 
that the paths which have the greatest survival probabil- 
ity are those with the smallest width (defined as the dis- 
tance between the extrema of the path), i.e. the sequences 



7 



— ^^---^ 



t=4 

t=8 

t=12 

t=16 

t=20 

t=24 

t=28 



FIG. 2: Survival probability P{S,t\x) given that the particle 
is at site x at early times and with pL = pR = 0.5. 
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FIG. 3: Single particle survival probability and bounds with 
PL = Pr = 0.5. 



x{t) = (0,1,0,1,0,...) and x{t) = (0, -1, 0, -1, 0, . . J 



This result gives support to the supposition in section [II 
that staying still (i.e. a diflFusion constant D' = 0) gives 
rise to the greatest chance of survival. We also estab- 
lished this to be case for two-dimensional walks up to a 
time t — 12. 

It is a simple matter to use the algorithm presented 
above to find the probability P{x,t\S) for the particle 
to be at coordinate x after time t given that it has sur- 
vived. Then, an application of Bayes' theorem yields the 
more telling quantity P{S,t\x), i.e. the probability that 
the particle has survived to time t given that it ends at 
coordinate x. The resulting data are plotted in Fig. ||and 
one sees quite clearly that that the particle is most likely 
to survive if it is at the origin, at least for times t < 28. 
This figure provides further weight to our assertion that 
staying still is the best particle survival strategy. 

As stated in the previous section, one can obtain esti- 
mates of various quantities at later times if one performs 
Monte Carlo sampling over particle paths. In fact, we 
produced histograms of P{S,t\x) this way and obtained 
data very similar to those shown in Fig. (except with 
poorer statistics). Hence we do not present them here. 
Instead we concentrate on the survival probabilities for a 
range of trap densities to compare with the bounds given 
by (H). 

First we consider the case of equal trap densities either 
side of the origin and the case of n = 1 and 2 particles 
starting at the origin. We generated the data for the 
case n — 1 using the algorithm described above, and 
densities pL — Pr = 0-5 until a time t = 30000. Bear- 
ing in mind the form of the bounds ( p2[ ) it is appropriate 

to plot the quantity X{t) ~ — h\Q(t) / ^ p^_Dt against log 

time. In all the simulations, D = D' = ^ and in this case, 
P- = 0.5. Hence the upper and lower bounds in ( |2^ ) con- 
verge to the constant A(oo) = 'i/y/n. Fig. || shows that, 
after an initial transient, A(t) does fall within the bounds. 
However, even at the late times probed in the simulation, 



— Numerical data 

— Bounds 



ln(p Dt) 

FIG. 4: Two particle survival probability taken from and 
bounds with pl ~ PR ~ 0.25. 



X(t) still seems to be far away from its asymptote. This 
highlights the fact that the predicted asymptotic form for 
the particle's survival probability (|l|) has not yet been ob- 
served in simulation, even with sophisticated methods at 
our disposal. 

The data for the case n = 2 have been taken from |l^ 
and are plotted with our bounds in Fig. ^. As with the 
case rt = 1 we have from ( p2[ ) that A(oo) — 4/^/?? and 
again the convergence to asymptopia is very slow. 

In Fig. 1^ we plot the single particle survival probability 
for the case where the densities of traps either side of 
the origin are unequal. Specifically we have the cases 
p+/p- = 2,4,8 with PR = 0.5 in each case. Note that 
the density used to scale the plots is always the smaller 
of the two, p- . Again we see that the numerical data lie 
within the bounds predicted by Eq. (p^. In these cases, 
however, the bounds we have presented do not converge 
so we have no predictions for the limiting value of A(t). 
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FIG. 5: Single particle survival probability and bounds with 
p+/p- —2,4: and 8. The symbols on the solid lines (represent- 
ing the numerical data) are included purely for the purpose of 
identifying each curve with the corresponding density ratio. 



V. ANALYTICAL RESULTS IN GREATER 
THAN ONE DIMENSION 

The upper and lower bounds on Q(t) derived in d = 1 
will be now be generalized to all d in the range 1 < d < 2 
and to c? = 2, the latter case requiring a slightly different 
treatment. The case d > 2 will also be discussed. 



A. Upper Bound 

Let the particle, with diffusion constant D', start at 
the origin, and the traps, with diffusion constant D, be 
randomly distributed in space with density p. As before, 
we assert, on intuitive grounds, that the "best strategy" 
for the particle is to stay at rest at the origin. With this 
assumption (which was verified numerically for d ~ 2, 
for all times up to t = 12, in the preceding section) the 
survival probability for D' = Q provides an upper bound 
on the survival probability for any D' > 0. Let Qi{t\r) 
be the probability that a given trap, starting a distance 
r from the origin, has not yet visited the origin at time 
t. It obeys the backward Fokker-Planck equation 



dQi 
dt 



D 



d~l dQi 



dr 



dr 



(33) 



where we have exploited the spherical symmetry of the 
problem. The boundary conditions arc Qi{t\0) = for 
all t and Qi{t\oo) = 1 for all t, while the initial condition 
is Qi{0\r) = 1 for all r > 0. Since there is no length scale 
in the problem, Qi{t\r) must have the scaling form 



Qiit\r) = fir/Vm) . 



(34) 



Substituting this form into Eq. ( ^3| ) gives an ordinary 
differential equation for f{x): 



dV d-ldf xdf 

— - H H = 

da;-^ a; da; 2 da; 



(35) 



with boundary conditions /(O) — 0, /(oo) = 1. The 
solution is 



2-d 



1 -1 /-a;"/* 



(36) 



For d = 1 our previous result, f{x) — erf(a;/2), is recov- 
ered. Note that Eq. ( ^6| ) is only valid for d < 2, since the 
integral diverges for d > 2. This regime will therefore 
require a different treatment. 

Eq. (|3^) gives the survival probability of a stationary 
particle in the presence of a single diffusing trap. Con- 
sider N traps in a large sphere of volume V centered on 
the origin. Each trap starts anywhere in the volume with 
equal probability. The average, over the initial positions 
of the traps, of the probability that none of the traps has 
yet reached the origin at time t is 



Q{t) = 



Vjv 



Dt 

1 - 



N 



N 



(37) 



Taking the limit N 
fixed, gives 



oo, y ^ oo, with p = N/V held 



Q{t) cxp 



1-/ 



(38) 



where the integral is now over all space. Inserting the 
function f{x) from Eq. ( ^6| ) and evaluating the integral 
gives the final result, which serves as an upper bound, 
Qu{t)i for the problem with general D' > 0: 



where 



Quit) = e^p[-adPiDtf^^] 



nd \ 2 



(39) 



(40) 



B. Lower Bound 

Our strategy for constructing a rigorous lower bound 
follows that employed in one dimension. We construct an 
imaginary (d-dimensional) sphere of radius I centered on 
the origin, and calculate the probability that (i) there are 
no traps inside the sphere at t = (ii) the particle stays 
inside the sphere up to time t, and (iii) no traps enter 
the sphere up to time t. As before, the set of trajectories 
(of particle and traps) selected by these constraints are 
a subset of all trajectories in which no traps meet the 
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particle, so the probability weight of this subset provides 
a lower bound on Q{t). We compute these probabilities 
in turn. 

(i) The probability that the sphere initially contains no 
traps is simply exp{- pYJ'') , where Vd = 271^^/^ /dr{d/2) 
is the volume of a d-dimensional unit sphere. 

(ii) The probability, Qp{t\r,l) that the particle stays 
inside the sphere up to time t is obtained by solving the 
backward Fokker-Planck Eq. (^3|), with D replaced by 
D', subject to the boundary conditions Qp{t\l, I) — and 
Qp(t\r,l) is analytic at r = 0, and the initial condition 
Qp{0\r,l) = 1 for r < I. The solution has the form 

oo 

Qp{t\r, l) = r''Y, eM-D'klt)J-,{knr) , (41) 

71=1 

where J-^{z) is a Bessel function of the first kind, 

,y = {2- d)/2 , (42) 

and fc„/ = Zn is the nth zero of J-^{z). The coefficients 
Cn are obtained from the initial condition, but their pre- 
cise values are of no interest here. Since the particle 
starts at r = 0, we need Qp{t\0,l). Its asymptotic form 



Qp{t\0,l) - exp{~^zfD't/f 



(43) 



(iii) To compute the probability, Qrit), that no trap 
enters the sphere up to time t (the target annihila- 
tion problem) we begin by calculating this probability, 
Qi{t\r,l), for a single trap. Then the probability that 
none of the traps enter the sphere is given by a natural 
generalization of Eq. (p§) , 



Qrit) = exp 



f d''r{l-Qi{t\r,l)} 

Jr>l 



(44) 



In contrast to the case I — used for the upper bound, 
there is no simple scaling form analogous to (|3^ ) for 
Qi{t\r,l) because I provides an additional length scale. 

The function Qi{i]r,l) obeys the backward Fokker- 
Plank equation (|33|), with boundary conditions 
Qi{t\l,l) = for aU t, Qi{t\(X),l) ^ 1 for aU t, 
and initial condition Qi{{)\r,l) = 1 for r > I. The 
solution can be found by Laplace transform techniques. 
The resuh is isl 



9 /r\v r°° elk 
(ik,/) = -(y) —exp{-DkH)G,{kr,kl), 



where 



(45) 



and Y^{z) is a Bessel functions of the second kind. 

Before continuing, we can first simplify Eq. (jij) as 
follows. First define 



r>l 



d'^r {1-Q,{t\rj)} , 



(47) 



where F{0) = follows from the initial condition 
Qi(0|r, Z) = 1 for all r > I. Then we use the backward 
Fokker-Planck equation (pSj) to write 



dtF 



-D 



d''rW^Qi{t\r,l) 



r>l 



= -D [ dA-S/Qi{t\l,l) 

J A 



^ DSdl^-^drQl{t\r,l)\r=l, 



(48) 



where A is the surface of the sphere, dA is a surface 
element directed along the inward normal to the sphere, 
and 



Sd — 



27^ 

T{d/2) 



(49) 



is the surface area of the unit sphere in d dimensions. In- 
tegrating the result (^) with respect to time, with initial 
condition F{0) = 0, Eq. (ph takes the form [El 



]T{t) = exp 



-pDSdr 



dt'drQlit'\r,l)\r=l 



(50) 

We are interested in the behavior of Qrit) for large t. At 
this point it is convenient to discuss separately the cases 
1 < d < 2, d = 2, and d > 2. 



1. The case 1 < d <2 

For 1 < d < 2, the hmction Qi{t\r,l), given by Eq. 
(Eq), has the large-t expansion |23] 



Qi{t\r,l) = 



- 1 



r(i + i 



2v) 



(51) 



r2(i + j.)r(i 

where t = ADt/P and we recall that — {2 — d)/2. 
Taking the derivative with respect to r, setting r = I, 
inserting the result into Eq. (pO[), and evaluating the in- 
tegrals over t' , gives the probability, that the target has 
not been annihilated by a trap, 

Qrit) = exp[^adp{Dtf/^-bdp{Dtr 

where ad is given by Eq. ( ^0| ) and 

_ 2^'^-^TT^/^r{d/2) 



\d-lp-d 



(52) 



(2-d)r2(i-d/2)r(d) 



(53) 



Note that, as with the case d = 1, the leading term is 
independent of I. This phenomenon can be attributed to 
the recurrence of diffusion in dimensions d < 2. 

Finally we assemble the contributions (i)-(iii) above to 
obtain a rigorous lower bound on the asymptotic behav- 
ior for 1 < d < 2, 



Xt) ^ eM-adp{DtY'^] X eM-pVdl" 



(54) 
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As usual, for a given time t* we choose a sphere radius I* 
to optimize the lower bound. The dominant ^-dependent 
terms for t ^ oo are the final two terms in the second 
exponential. Ignoring constants of order unity, we find 
that the value of I* that gives the greatest lower bound 
is 

l/(4-d) 

I* - ( 1 (i*)(2-rf)/(4-rf) (55) 



Inserting this into (^4|) the second exponential takes the 
form 

exp[-const(i?')''"''/''~'^(p^'"')'^^'"'^^'/''"'^]- (56) 

The neglected first term in the second exponential in ( |5^ ) 
behaves as l'^ ~ t'^(2-t')/(4-<i)^ which is indeed negligible 
compared to for large t (recalling that d > 1 

here). 

In summary, the best lower bound behaves as 

QL{t) ^ exp[~adpiDtf/^ + 0(1"/^^-^^)]. (57) 

Since (i/(4 — d) < d/2 for d < 2, the two bounds pinch 
asymptotically, to give the exact result 

where we recall that ad is given be Eq. (^ 
Ad in Eq. (§) is therefore given by 



The constant 



A. 



l<d<2 



(59) 



Note that the subdominant term in ( pT\ ) decays more 
slowly relative to the leading term as d ^ 2, signaling a 
change of behavior at c? = 2. Note also that the coefficient 
ad vanishes at d = 2, suggesting a slower decay than a 
simple exponential in two dimensions. We now show that 
this expectation is correct, and determine the constant 
A2 in Eq. d). 



where zi is now the smallest zero of Jo{z). The dominant 
terms in the exponential for large t are the first and last 
terms. Extremizing this bound with respect to / at some 
fixed t* gives 



to leading order, and 



Ql(^) ~ cxp 



AnpDt 



\n{pDH/D') 



O 



tln(lnt) 
In^t 



(63) 



(64) 



As far as the upper bound is concerned, Eq. ( p9[ ) is not 
useful in d = 2, since 02 = 0. This tells us that the prob- 
ability that a trap will reach a specified region of zero 
volume (i.e. a specified point) is zero in two dimensions. 
The A particle has to be given a non-zero size (or the 
system put on a lattice) for a non-zero trapping proba- 
bility. We therefore assign the particle a non-zero radius 
a, but still treat it as stationary for the upper bound. 
The traps will, for the moment, continue to be treated as 
point particles. With the definition that trapping occurs 
if a trap enters within the particle's radius (so that a is 
an interaction range), our upper bound is just given by 
the probability Qrit), Eq. (61), but with I replaced by 
a, to give 



Quit) ~ exp 



AnpDt 



O 



(65) 



In the limit t — > 00, the bounds converge to give the 
asymptotic result Q{t) ^ exp(~4:Tr pDt / Int) or, equiva- 
lently 



IntlnQ(i) ^ 
hm = 47r . 

*— >oo pDt 

This gives the constant A2 in Eq. (M as 



(66) 



2. The case d=2 



For d — 2 the asymptotic form of Qi(t\r, I) is pi 



Qi{t\r, I) = 2\n{ J 



1 

Inr 



O 



1 



In^ T 



(60) 



where t = ADt/P as before. Inserting this into Eq. (50), 
with d = 2, gives the probability that no trap has entered 
the circle of radius I up to time t: 



?t(0 = exp 



AnpDt 



O 



t 



(61) 



ln(4i:ii/'^) 

Following our previous procedure, the asymptotic lower 
bound is given by 

A-KpDt ^ ( t 



At) 



exp 



ln(4Dt/?2) 
,2 2D't 



O 



(62) 



A2 = AnpD 



(67) 



As noted previously, the algorithm described in Sec- 
tion 1^ can be used to simulate the trapping reaction 
in any integer dimension. Numerical results for the two- 
dimensional system were presented in 16 and we com- 
pare these data with the asymptotic result ( p6[ ) in Fig. ||. 
Wc find that the deviation of the numerical results from 
the asymptote is even more marked in two dimensions 
than in one (see Fig. Part of the reason for this 
is, presumably, that the increased number of sites in 
two dimensions means that one cannot probe such late 
times. A second, and perhaps more important, reason 
is the very large corrections to scaling evident from our 
bounding arguments. The relative size of the sublead- 
ing term in Eq. ( |65| ) is 0(1/ hit) while the subleading 
term for the lower bound, Eq. (|64|), is even larger, of 
relative size 0(ln[lnt]/ Ini). This suggests that conver- 
gence to asymptopia will be extremely slow in two dimen- 
sions. Note also that the particle's survival probability 
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FIG. 6: Numerical data for the two-dimensional trapping re- 
action taken from In the simulation, the trap density 
p — 1/4: and the diffusion constants D — D' = 1/4. The 
asymptote given by Eq. (pfj) is plotted for comparison. 



was found to decay to ~ 10~^^ after t = 1600 time steps. 
This emphasizes the importance of determining the cor- 
rections to asymptopia in order to determine the form 
of the survival probabihty in numerically (and, indeed, 
experimentally) accessible regimes. 



3. The case d > 2 

The same bounding arguments can be applied equally 
well in d > 2. The main difference from d < 2 is that the 
bounds no longer converge, so it is not possible to deter- 
mine Xd exactly (except for d very close to 2 — see below) . 
The basic idea is the same as for d < 2, except that the 
particles must be given non-zero sizes (or, equivalently, a 
non-zero range of interaction). We let the particle have 
radius a, and the traps radius b. A reaction is deemed to 
have occurred if there is an overlap between the particle 
and any trap, i.e. if the centers approach more closely 
than a distance R = a + h, which is the range of inter- 
action. (Note, however, that we continue to assume that 
the traps do not interact with each other. In particular, 
there is no excluded volume interaction between traps.) 

The upper bound is obtained from the target annihila- 
tion problem with target radius R. For d > 2, the single- 
trap survival probability, Qi{t\r, R), has a non- vanishing 
large-t limit given by the well-known result 



Qi{oo\r,R) = l-{R/r) 



d-2 



(68) 



which is easily obtained from Eq. (33) on setting the 
left side to zero, and imposing the boundary conditions 
Qi{oo\R,R) = 0, (3i(oo|oo, i?) = 1 on the resulting or- 
dinary differential equation. Inserting this form in Eq. 
(pO|), with / = R, and evaluating the time integral, gives 
an upper bound with the leading large-t behavior 



Quit) exp[-(d - 2)SdpR''-^Dt] 



The lower bound is obtained in a similar fashion, fol- 
lowing the pattern established for d < 2. One constructs 
a notional sphere of radius Z, centered on the initial posi- 
tion of the particle. The bound is given by the subset of 
trajectories in which (i) there are no traps initially within 
the sphere, (ii) the center of the particle remains within 
a sphere of radius / — a, so that the particle remains en- 
tirely inside the sphere of radius I, and (iii) the center 
of every trap remains outside a sphere of radius I + b, 
so that every trap remains entirely outside the sphere of 
radius I. 

The probability of (i) is eyi^[-pVd{l + hY]. The proba- 
bility of (ii) has the asymptotic form exp[— zJZ3't/(Z— a)^], 
where zi is the first zero of J-y{z) = J(d-2)/2{^- The 
probability of (iii) is given, for large i, by Eq. (|69|) with 
R replaced by / -I- 6. Assembling these three contributions 
gives the asymptotic lower bound 



QL{t) ^ exp[-(d-2)5dp(r 



bf-'^Dt- zlD't/{l-af 
~pVd{l + bY] . (70) 



This has to be maximized with respect to I. For i oo, 
the first two terms in the exponent dominate, and the 
final term is negligible. Setting / = a+x, and maximizing 
with respect to x, gives the equation 

(d - 2fSdpD{x + Rf-^ = 2zfD'/x^ , (71) 

where R = a + b as before. This equation cannot be 
solved analytically for general d, so we concentrate on 
two soluble cases — the physically interesting case d = 3, 
and the limit d — > 2+. 

For d = 3, we have Sd — 47r and zi — tt, giving x — 
{TTD'/2pDy/^ and 



Ql(0 - exp[-4:TTpDRt - 'd{2Ti'^ pD^ D'f'^t 



(72) 



Combining the two bounds, we obtain the asymptotic 
form Q{t) ~ exp(— Aat), as in Eq. (P, with the bounds 

AirpDR < A3 < AupDR + i{2Ti^ pD^/lyfl^ , d = 3 . 

It is worth noting that the second term on the right is 
negligible compared to the first if D'/D <C pR^ , i.e. when 
D' / D is small compared to the number of traps per in- 
teraction volume. 

For d = 2 -t- e, Eq. ([7l|) has the solution 



1 (4D'\ 
e V T^pD J 



1/2 



(74) 



to leading order for e ^ 0, giving the lower bound 
Qbit) ^ exp[— 27rep£'i] to leading order in e. In the same 
limit, the upper bound (^9|) has exactly the same form, 
giving the result Q{t) ^ exp(— A^i) with 



Xd = 2iTpD e 



(75) 



(69) 



Hence the bounds pinch to leading order in e, but not for 
general d. 



12 



To conclude this section, we consider again the case 
where n particles start from the origin, and we want the 
probability that all survive until time t. As noted in 
the discussion of the one-dimensional case, n only enters 
in the calculation of the lower bound, in the term giving 
the probability for the particle to stay inside the notional 
box (d = 1), or sphere (d > 1), of size This probabil- 
ity behaves as exp(— const£>'t//^), so having n particles 
simply requires raising this factor to the power n, which 
is equivalent to replacing D' by nD' . Since the asymp- 
totic forms we derive do not depend on D' for d < 2, 
it follows that our results are independent of n in this 
regime. For d > 2, however, our results do depend on D' 
(see Eq. ([73|)). In this regime, therefore, the generaliza- 
tion to arbitrary n is achieved through the replacement 
D' nD'. 



VI. DISCUSSION AND SUMMARY 

In this paper we have derived a number of results for 
the asymptotic survival probability of a particle diffusing 
among randomly distributed diffusing traps with density 
p. We allow the particle and traps to have different diffu- 
sion constants, D' and D respectively. Our results take 
the forms originally derived by Bramson and Lebowitz 
[0, as expressed in Eq. With one assumption, sup- 
ported by numerical evidence, we have obtained exact 
results for the coefficients Xd in (|l]) for dimensions d < 2, 
and an exact inequality for dimensions d > 2. These re- 
sults are given by Eqs. (|l|), (H), and (|7|). For d < 2 the 
results for are independent of the diffusion constant 
D' of the particle. 

The results are obtained by deriving upper and lower 
bounds for A^, and showing these coincide for d < 2. 
Whilst our lower bound is rigorous, we had to assume 
that the particle's survival probability for D' = pro- 
vides an upper bound on its survival probability when 
D' > when the trap distribution is symmetric. Indeed, 
for the d — 1 system with different densities of traps to 
the left and right of the particle, it was found that staying 
still is not the particle's best strategy. Instead, trajecto- 
ries that survive for long times tend to be those in which 
the particle drifts to the side with the lower trap density. 
This emphasizes the crucial role of the symmetry of the 
trap distribution, an observation supported by perturba- 
tive studies for a system with a finite number of traps 

In all cases the particle and traps are assumed to move 
in a continuous space, and to have zero size for d < 2. 
For d > 2 is is necessary for the particle and/or the traps 
to have non-zero size (otherwise the survival probability, 
for motion on a continuous space, is one for all time). We 
also take the traps to be randomly distributed in space at 
time t = 0, with uniform density p. This raises the ques- 
tion of the extent to which the results are "universal" , 
i.e. independent of the microscopic details of the model, 
a question which we now address. 



We argue that, for d < 2, the results do indeed have 
a degree of universality. In d = 1, the optimal box size 
used to obtain the lower bound on Q{t) is large, I ~ t^^^, 
as t — *■ oo, so the effect of the particle having a finite 
size when confined to this box is negligible. With a little 
thought one sees that the same is true for all d < 2. The 
optimal length scale for the lower bound grows with time 
as I ~ i(2-d)/(4-d) (d < 2) or Z ~ Ini (d=2), and the 
results are independent of the particle and trap sizes, as 
far as the leading-order results are concerned. The same 
is true of the upper bound — the finite-size corrections 
come in at subleading order. 

The dominance of large length scales at late times, for 
d < 2, also suggests that the asymptotic results are inde- 
pendent of whether the model is defined on the contin- 
uum (as here) or on a lattice, an assumption implicitly 
made earlier when we compared our theoretical predic- 
tions to numerical results obtained from lattice simula- 
tions. For d > 2, however the dominant value of I that 
determines the lower bound is time-independent. There- 
fore we expect a lack of universality in this case. The 
explicit dependence on the interaction range R in Eq. 

) is a signature of this effect. Note, however, that to 
leading order in e = d— 2, Ad is independent of R (see Eq. 

)) and we expect the result to be universal to this or- 
der. Physically, this is because the length scale I = a + x 
diverges as e ^ (see Eq. (]7^)). 

A further universality question concerns universal- 
ity with respect to the initial conditions. We have 
taken Poissonian initial conditions, where the probabil- 
ity P/v(y) of having N traps in a volume V is given by 
Pn{V) = [{pV)^ /N\]exp{^pV) for any V. The lattice 
simulations, where the number of traps on each site has a 
Poisson distribution (with mean p, say) has this property, 
namely the number of traps on m sites has a Poisson dis- 
tribution with mean mp. Whether there is a larger class 
of initial conditions sharing the same asymptotic behav- 
ior is a question deserving further study. 

We conclude by discussing some recent papers related 
to the present work, and directions for future work. The 
coefficient Ad in Eq. (P has recently been calculated using 
a diagrammatic method |^ to first order in (2 — d). The 
quoted result, however, exceeds our rigorous upper bound 
for Xd (corresponding to the lower bound for Q{t)) by a 
factor of two. This is because in |^ Ad depends on D 
and D' only through their sum D + D' psf whereas our 
rigorous upper bound on Ad depends solely on D. It is 
interesting to note that in the related process A + B ^ 
A where the single A particle acts as a trap for the B 
particles, certain properties of the B-particle distribution 
can be expressed as functions of D + D' However, 
we stress that for the A + B ^ B reaction studied here, 
the asymptotics are entirely governed by the i?-particle 
diffusion constant for d < 2. 

In a very recent work [^0| our approach, as outlined in 
p5[ , has been generalized to diffusion on fractals for the 
case where the fractal dimension of the traps' trajectories 
is greater that the physical dimension (this condition is 
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the analog of the condition d < 2 in the present work). 
It should be noted, however, that in Ref. |Q the opti- 
mal lower bound on Q{t) is not obtained. For d < 2, 
only the subdominant corrections to the leading terms 
are affected, and the upper and lower bounds still pinch 
asymptotically. For d = 2, however, the approach used 
in yields bounds that no longer converge at large i, 
so the exact result ( |67| ) for A2 is missed. 

In this paper and our earlier work we noted that 
the extant simulation data (l6| fail to reach the asymp- 
totic regime even though survival probabilities are so 
small that they can only be measured using sophisticated 
methods. Perhaps the most important challenge, there- 
fore, is to obtain a better understanding of the correc- 



tions to asymptopia in order to make testable, quanti- 
tative predictions. Other directions for future work in- 
clude exploring further the extent to which our results 
are universal, and establishing rigorously the validity of 
our upper bound. 
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